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We report the first low temperature magnetotransport measurements on electrochemically fab- 
ricated atomic scale gold nanojunctions. As T ^ 0, the junctions exhibit nonperturbatively large 
zero bias anomalies (ZBAs) in their differential conductance. We consider several explanations and 
find that the ZBAs are consistent with a reduced local density of states (LDOS) in the disordered 
metal. We suggest that this is a result of Coulomb interactions in a granular metal with moderate 
intergrain coupling. Magnetoconductance of atomic scale junctions also differs significantly from 
that of less geometrically constrained devices, and supports this explanation. 

PACS numbers: 



At low temperatures disordered metals exhibit a re- 
duced local density of states (LDOS) at the Fermi 
level, seen as a zero bias anomaly (ZBA) in tun- 
neling spectroscopy {e.g. 0). This is the result 
of disorder-enhanced electron-electron interactions 0- 
Granular metals with large intergrain conductances {g = 
(Gjj)/(2e2//i) >> 1) are predicted 3J to have LDOS 
suppressions approaching the perturbative result Q for 
weakly disordered films. With strong disorder and geo- 
metric constraint j£ i5j I the ZBA in metal films can ap- 
proach 100%, ascribed to a correlation gap due to strong 
Coulomb interactions^. Similarly, weakly coupled gran- 
ular metals {g « 1) should act as arrays of tunnel 
junctions, with an exponentially suppressed tunneling 
LDOSfl as T ^ 0. 

Metallic nanojunctions (MNJs) are tools to exam- 
ine geometrically constrained, disordered metals on the 
nanometer scale. While clean break junctions made in ul- 
trahigh vacuum (UHV) have been studied extensively 0, 
nanojunctions between highly disordered metals are com- 
paratively unexplored. 

We present the first low temperature measurements 
of atomic scale metal junctions made by electrochemical 
deposition, a method proposed for molecular electron- 
ics investigations 0, El- Such metals may be disordered 
by grain boundaries, ionic impurities, and surface adsor- 
bates. As we reported elsewhere, lOj, in "large" junctions 
(G(300 K) >> Go = 2eV/i), small ZBAs are consis- 
tent with the perturbative theory of Altshuler, Aronov 
and Lee (AAL)[2. When G(300 K) - 2e^/h, how- 
ever, low T conductance suppression approaches 100%. 
Here we consider several models and show that these 
junctions are atomic scale probes of the LDOS of the 
disordered metal leads, which exhibit non-perturbative, 
temperature-dependent LDOS corrections. We give a 
phenomenological description of the LDOS suppression, 
and argue that its likely origin is the granular charac- 
ter of the electrodeposited material. Finally, these junc- 
tions exhibit nontrivial weak localization magnetocon- 
ductances, consistent with junction size and the hypoth- 
esis of granularity. 

The MNJs are prepared by electrochemistry starting 
from electrodes defined by e-beam lithography. The de- 



FIG. 1: Schematic of the electrode configurations used to 
make the nanoscale junctions. At right, overall electrode ge- 
ometry; at upper left, pre-deposition image of tip geometry 
used for sample C; lower left, pre-deposition image of tip ge- 
ometry used for samples A and B. Samples A, B, and C are 
representative atomic scale junctions, as discussed below. 

tails of sample fabrication are available elsewhere p^. 
Gold electrodes 25 nm thick with separations of ~20 nm 
are prepared on p+ silicon substrates coated with 200 nm 
of thermal oxide. Electrode configurations are shown in 
Fig. ^ Lumped capacitance between each electrode and 
the substrate is estimated to be ~ 50 pF. The electrodes 
are covered by 20 nm AI2O3 during evaporation, limiting 
electrochemistry to the electrode edges. The evaporated 
Au has a typical resistivity of 5 /iil-cm. 

Additional gold is deposited using a buffered aqueous 
solution of potassium cyanaurate'sj, while interelectrode 
conductance is monitored with standard lock-in tech- 
niques. Discrete conductance steps on the order of Go 
are observed during junction formation, corresponding to 
atomic reconfigurations. The MNJ is grown to a specified 
conductance, rinsed in deionized water, and dried with 
dry nitrogen. Through measurements on test structures, 
we find that the average resistivity of the electrodeposited 
gold is ^ 35 iiri-cm at 4.2 K, corresponding to an elastic 
mean free path of ~ 2.5 nm, much shorter than that of 
the evaporated Au. 

Spontaneous conductance switching and its strong sup- 
pression as T is lowered suggest that the junctions con- 
sist of a small number of atoms that can diffuse read- 
ily at room temperature. Surviving junctions with G ~ 
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FIG. 2: Zero bias resistance as a function of temperature 
for three samples with room temperature conductances on 
the order of Go- The form of the dependence is nonac- 
tivated. Solid lines are fits: (A) 0.61 exp(1.91T"°-^); (B) 
0.91exp(8.06r-° **); (C) 1.65 exp(3.6ir-"-^^). 

1 — lOGo are stable for tens of minutes, typical of junc- 
tions not prepared in UHV 11]. 

We have successfully measured ten nanojunctions with 
room temperature conductances ranging from 0.5 to 200 
Go in a variable temperature cryostat. Using standard 
quasi-4-terminal lock-in techniques we have measured 
nanojunction differential conductance G{V,T) — dl/dV 
(and differental resistance) as a function of temperature, 
dc bias voltage, gate voltage, and magnetic field. All 
nanojunctions are Ohmic up to 200 mV at 300 K. 

As discussed elsewhere Il0|, high conductance junctions 
{e.g. G(300 K)- 30 Go)' typically exhibit a smaU (15% 
at 1.8 K) ZBA, logarithmic in temperature below 30 K. 
We interpret this ZBA and its scaling with bias voltage 
as consistent with the AAL perturbative DOS suppres- 
sion. An analysis of the Coulomb interaction with a sin- 
gle coherent scatterer[13j can quantitatively explain this 
data0|, provided the disordered metal is granular, with 
an intergrain conductance on the order of 50 Go; such an 
analysis is only applicable for systems with G >> Go. 

In atomic scale junctions prepared as above with 
G(300 K) ~ IGo, we find ZBAs approaching 100% con- 
ductance suppression as T ^ 0. Such deep ZBAs are 
never seen in the larger junctions and are present in both 
differential conductance and resistance measurements. 
The only difference between large and atomic scale junc- 
tions is a brief amount of electrodeposition time; the mi- 
crostructure of the deposited material is identical. Figure 
|2]shows the zero bias resistance vs. T on Arrhenius-style 
axes for three of these samples. 

We now present evidence that these large ZBAs are 
caused by nonperturbative low temperature corrections 
to the LDOS of the electrodeposited material. Fig- 
ure El (left inset) shows G{V — 0,T) for sample A 
(G(T=300K)=lGo) as a function of temperature cy- 
cling. During initial cooling (trace 1), G varies little 
until T < 50 K, when the zero bias suppression begins. 
At 15 K, G spontaneously increases by about 0.6 Go. 
The conductance then continues to decrease upon cool- 
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FIG. 3: Zero bias anomaly for Sample A at (top to bot- 
tom) 20 K, 10 K, and 2 K. Solid lines are the phenomeno- 
logical model described in text. Left inset: various branches 
of G(0, T) for this sample, each branch corresponding to a 
different arrangement of contact atoms. Right inset: all the 
branches collapse onto a single curve when each is scaled by 
a multiplicative constant. 

ing. When the nanojunction is warmed (trace 2), the 
high temperature conductance appears to be ~ 2Go. At 
220 K, however, G spontaneously decreases by about 0.85 
Go, returning near its original value. Repeated ther- 
mal cycling and LED illumination at 2 K result in other 
branches of G vs. T (traces 3,4,5). these G(y = 0, T) 
curves collapse onto one curve (right inset) when each 
branch of G{T) is multiplied by a non-integer constant. 

The discrete changes in G strongly support the idea 
that this junction's low room temperature conductance 
(^ Go) is due to the junction's atomic scale. The addi- 
tion of a single partially transmitting channel upon cool- 
ing occurs as thermal contraction slightly decreases the 
interelectrode distance. Thermal expansion on warming 
stretches the junction, and the additional channel is lost, 
just as in hysteresis seen in mechanical break junction 
measurements^^. These serial rearrangements indicate 
that this junction was in a fortuitous regime of stability. 
In the other samples, when rearrangements took place 
the junctions either broke completely or coalesced to a 
high conductance state. 

For each conductance branch, the ZBA is measured at 
several temperatures. Figure |21 shows the conductance 
versus bias voltage at 2 K, 10 K, and 20 K for this nano- 
junction in one of its configurations. The same factors 
used to collapse the G{V = 0, T) branches also collapse 
the bias sweep data onto a single set of curves. It is clear 
that G(y, T) is only multiplicatively scaled by discrete 
atomic rearrangements of the small number of conduct- 
ing channels. 

We now consider possible origins of such a nonpertur- 
bative ZBA. A successful explanation must be consistent 
with: (a) the temperature dependence of the zero bias 
conductance; (b) the functional form of the ZBA vs. bias 
voltage and temperature; and (c) the scaling of the ZBA 
data with junction conductance as shown in Fig. |3| 

Coulomb blockade in an array of grains weakly coupled 
by tunnel junctions can produce a ZBA 16]. However, the 
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FIG. 4: Bias sweep data for two other low conductance junc- 
tions. Temperatures for sample B (left) are, top to bottom, 
40 K, 20 K, 10 K, 4 K, and 2 K. Temperatures for sample C 
(right) are, top to bottom, 20 K, 10 K, and 2 K. Solid lines 
are calculated using the phenomenological LDOS described in 
the text. 

temperature dependence in Fig. 13 and the ZBA width 
of our data are incompatible with such a simple model. 
For a small number of metal islands, one would also ex- 
pect asymmetries or Coulomb staircase effects, neither of 
which are observed. Gate voltage sweeps over a 100 V 
range produce no measurable response at any tempera- 
ture considered. Finally, this picture would not explain 
the scaling in Fig. |31 

Another candidate is environmental Coulomb blockade 
(ECB)^3|- Conduction through a single junction in se- 
ries with an environmental impedance is suppressed at 
low T and V because tunneling electrons excite electro- 
magnetic modes of the environment, reducing the phase 
space for tunneling^^. The IV characteristics are deter- 
mined by P{E), the probability of a tunneling electron 
to excite an environmental mode of energy E. The form 
of P{E) depends on the impedance of the environment, 
which typically must be well controlled to observe ECB 
in single junctions. 

ECB has been seen^^ in clean point contacts made 
using mechanical break junctions. The environmental 
impedance was modeled as the junction capacitance in 
parallel with the lead resistancefigj. Our voltage and 
temperature scales would imply a charging energy 10^ x 
that in Ref. |l9l |. This seems unlikely given our similar 
geometry and greater stray capacitance. Furthermore, in 
atomic scale junctions the ECB suppression scales with a 
Fano factor under changes of channel transmittance j2(il| . 
contradicting Fig. |3 This Fano scaling is expected to 
hold for a general environmental impedance 2 1], seem- 
ingly ruling out ECB as the origin of the observed ZBAs. 

A tunneling model'5| is a natural approach to ana- 
lyzing the ZBA data. One can consider a single tunnel 
junction with a transmission probability T(E) and an ef- 
fective LDOS i'{E, T). The conductance G normalized by 
the high bias / high temperature conductance Gh is then 
proportional to T{E)v{E,T)/vq convolved with a ther- 
mal spreading function. Here is the background effec- 
tive density of states. Our ZBA data exhibit a consider- 
ably stronger T-dependence than the spreading function; 
therefore a T-dependent LDOS is required in this model. 



In this picture the transmittance is independent of T and 
(slowly varying in) E, but depends on the precise config- 
uration of the few atoms that make up the junction. The 
LDOS is determined by disorder "built in" during the 
electrodeposition process. This would naturally explain 
the scaling in Fig.|3 different branches are equivalent to 
probing the same LDOS suppression using tunnel junc- 
tions of various sizes. 

The large ZBAs in the atomic scale junctions, rem- 
iniscent of correlation gaps in highly disordered metal 
films|3,l3, imply a nonperturbative LDOS suppression. 
Efetov and TscherichJ3| consider the LDOS in granular 
metals having a dimensionless intergrain conductance g. 
They derive expressions for the LDOS in the limit g » 1 
{g << 1), but the calculated T-dependence is too gradual 
(steep) to fit our data. This suggests that our samples 
fall between these extremes, into the intermediate range 
of g for which there is currently no analytic expression. 

A generalized treatment of corrections to electron 
tunneling|22j| accounts nonperturbatively for both intra- 
electrode Coulomb effects (the AAL LDOS correction) 
and interelectrode Coulomb interactions in the presence 
of an electromagnetic environment {P{E) theory). Ef- 
fective tunneling densities of states have been calculated 
that agree well with experiments on spatially extended 
tunnel junctions '2^. One would expect some form of 
this generalized theory to apply in the granular metal 
case. We introduce an ansatz for the functional form of 
5i^(e,r): 

Si^je, r) ^ zyp I 1 - erf, / — '""^ | (1) 

where F, a, and b are sample-specific, temperature in- 
dependent parameters. This model LDOS is able to de- 
scribe empirically the ZBA data over a broad tempera- 
ture and voltage range. The T = version of this form 
is derivedp^ for Id tunnel junctions, and is also equiv- 
alent to ECB in an ultrasmall junction connected to an 
RC transmission line [l8Ll2^ . However, the Fano factors 
discussed above for ECB in point contacts make this in- 
terpretation difficult to reconcile with the observed scal- 
ing. The solid lines in Figs. O and ^ are fits using Eqn. 
(Q. The relevant parameters is shown in Table 



Sample 


Gh [2e^/h] 


er [J] 


a [J] 


b [J/K] 


A 
B 
C 


1.51* 
1.29 
0.57** 


2.96 X 10"^^ 
5.11 X 10"^^ 
6.07 X 10"^^ 


2.33 X 10"^^ 
1.70 X 10"^^ 
1.87 X 10"^^ 


1.43 X 10"^^ 
6.94 X 10"^^ 
1.57 X 10"^^ 



TABLE I: Model parameters from Eq. used to reproduce 
the ZBA data for the samples shown in Fig. |21 *For sample 
A, data from branch number 3 from Fig. El were used for the 
fits. **For sample C, a slight junction rearrangement led us 
to use Gh = 0.62 to fit data at 2 K, with the other parameters 
unchanged. 

We also observe magnetoconductance (MC) in atomic 
scale junctions that differs from that in larger, cleaner 
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FIG. 5: Weak localization magnetoconductance sweeps at 
several temperatures for sample B, the junction from Fig. 2] 
Curves have been offset vertically for clarity. 

gold systems. Quantum MC effects are dominated by the 
coherent volume centered on the conductance-limiting 
junction. Larger junctions and lead electrodes exhibit 
weak antilocalization similar to that observed in clean 
gold films, because the spin-orbit scattering time, rgo, is 
more rapid than the inelastic time r<^|23|. While rso is 
temperature independent, typically increases with de- 
creasing T, and > rso for evaporated Au films below 
20 K. 

Figure shows MC data for Sample B (data for other 
atomic scale samples are similar). Curves are offset by 
0.0005Go- Positive MC contributions suggest that 
is comparable with rgo, with apparently rgo > for 
T >- 6 K. 

Quantitative analysis is difficult because the coherent 



volume around the constriction is of uncertain dimen- 
sionality with respect to diffusion and quantum coher- 
ence. As a qualitative guide, solid curves are fits using 
Id weak localization, assuming 3d diffusion, an effective 
width of 3 nm (roughly sets the field scale) and an ef- 
fective length of 5 microns (essentially a fixed numerical 
factor to scale the conductance axis), and allowing 
and Lso to vary at each temperature. A natural ex- 
planation for a varying Lso would be a fixed rso but 
a temperature dependent diffusion constant, D{T), that 
decreases as T — > 0. Such a dependence supports the 
granular metal model, where D{T) would vary as inter- 
grain interactions become important with decreasing T. 
We find that the ratio (L^/Lso)^ = t^{T)/tso (indepen- 
dent of D) increases with decreasing T roughly like l/T. 
The fact that r^/rso ^ 1 here suggests that is shorter 
here than in larger junctions. 

These first low-T studies of electrochemically made 
atomic scale Au nanojunctions have revealed signifi- 
cant departures from the properties seen in larger and 
cleaner Au structures. Atomic scale junctions locally 
probe a suppressed density of states in the disordered 
leads. We provide an expression that describes the data 
pheonomenologically, and suggest an underlying physical 
origin for the effect consistent with magnetoconductance 
data. Further studies are required to determine conclu- 
sively the physics behind the gap observed in the LDOS 
and the nontrivial MC properties of this system when 
examined at the atomic scale. 
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